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O 

o, 

0^ , We study a model of small-amplitude traveling waves arising in a supercritical Hopf- 

bifurcation, that are coupled to a slowly varying, real field. The field is advected by 

f^ ' the waves and, in turn, affects their stability via a coupling to the growth rate. In the 

absence of dispersion we identify two distinct shortwave instabilities. One instability 

induces a phase slip of the waves and a corresponding reduction of the winding number, 

C ' while the other leads to a modulated wave structure. The bifurcation to modulated 

K^ . waves can be either forward or backward, in the latter case permitting the existence of 

►v> , localized, traveling pulses which are bistable with the basic, conductive state. 



1 Introduction 

Various systems have been shown to exhibit bifurcations to spatially periodic patterns in the 
presence of an additional, weakly damped mode. If that mode corresponds to a conserved 
quantity it constitutes a true zero mode as in, e.g., translation of the fluid interface in two- 
layer Couette flow |I[] or a shift in the displacement velocity of seismic waves in a viscoelastic 
medium . Examples of weakly damped modes include the large-scale concentration field in 
binary-fluid convection |^, ^ or the real, slow mode in the 4-species Oregonator model of the 
BZ reaction ||^. The evolution of instabilities in these and similar systems will be coupled 
to the dynamics of the weakly-damped mode. Within the context of a weakly nonlinear 
approach near onset of the pattern-forming instability, the Ginzburg-Landau model is altered 
through the coupling to an evolution equation for the additional mode. The precise form of 
the coupling depends on the physics and symmetries of the system. 

Steady instabilities in conserved systems with and without reflection symmetry have been 
considered. With reflection symmetry the instability, setting in at finite wavenumber, was 



found by Mathews and Cox to be amplitude-driven, leading to a supercritical modulation 
of the underlying pattern. In the absence of reflection symmetry, Malomed identified 
modulated patterns by deriving longwave equations for the phase of the underlying pattern 
coupled to the zero-mode. A particle-in-a-potential model in a moving frame was then 
derived for the local wavenumber, which admits both traveling and steady modulations. 
Ipsen and S0rensen ||^ investigated the effect of a real, slow mode in reaction- diffusion 
systems near a supercritical Hopf bifurcation. They found the slow mode leads to new finite- 
wavenumber instabilities which alter traditional Eckhaus and Benjamin- Feir stability criteria 
for periodic waves. Barthelet and Charru studied instabilities of interfacial waves in two- 
layer Couette-Poiseuille flow with Galilean invariance and no reflection symmetry, where 
the zero-mode is a shift of the fluid interface. The corresponding model of waves coupled 
to the zero- mode, derived and worked out by Renardy and Renardy [§], results in multiple, 
separated stability regions for periodic waves in contrast to the Eckhaus stable band exhibited 
by a single Ginzburg-Landau equation. The critically damped mode corresponding to a two- 
dimensional mean-flow in stress-free convection becomes relevant for not too large values of 
the Prandtl number. Bernoff [^ derived a Ginzburg-Landau equation coupled to the mean- 
flow mode from the Bousinesq equations and calculated the skew-varicose and oscillatory 
skew-varicose stability boundaries. This calculation brought into agreement previous work 
by both Siggia and Zippelius [^ and Busse and Bolton |Tl[] who investigated the effect of 
mean-flow modes with non-zero vertical vorticity on the stability of rolls in Rayleigh-Benard 
convection. 

In binary-fluid convection, Riecke [^ derived an amplitude equation for traveling waves 
arising in a Hopf bifurcation coupled to a critically damped mode related to large-scale 
modulations of the concentration field in the limit of small Lewis number. Instabilities of 
periodic waves within these equations were found to lead to localized traveling pulses which 
exhibited the anomalous slow drift observed in experiment |]T2[, and in numerical simula- 
tions |13|. Riecke and Granzow presented a similar model on phenomenological grounds to 
describe the dynamics of traveling waves in electroconvection in nematic liquid crystals [P^ . 
This model for oblique (zig and zag) waves coupled to a slow mode, possibly corresponding 
to a charge- carrier mode, exhibits instabilities to localized, worm-like structures as seen in 
experiment |]15|. A striking feature of the worms is that although the extended waves bifur- 



cate supercritically from the basic, conductive state, the worms themselves are bistable with 
the latter. 

In the present paper, we consider traveling waves arising in a supercritical Hopf-bifurcation 
coupled to a real, slowly varying field. The field is advected by the waves and, in turn, can 
affect the dynamics of the waves through coupling to their growth rate. We briefly introduce 
the model and present a cursory linear analysis in section 2 which reveals distinct phase and 
amplitude instabilities of the waves. In sections 3 and 4 we are concerned with characterizing 
the linear and nonlinear behaviors of the phase and amplitude instabilities, respectively. In 
both cases we derive an envelope equation for the instability near its threshold and compare 
the results to full numerical simulations of the original equations. The envelope equations, 
in agreement with the numerics, indicate the phase-instability leads to a backward Hopf 
bifurcation, while the amplitude-driven instability leads to modulated waves, arising super- 
or subcritically. In the latter case, the subcritical branch can be bistable with the basic, 
conductive state and localized wave pulses arise. 



2 The Extended Ginzburg-Landau Equations 

We consider a model of a traveling wave with complex amplitude A and group velocity s 
coupled to a real, weakly damped mode C The general form of the equations to orders 
e^, e^ respectively, where e is a measure of the distance from threshold of the bifurcation to 
traveling waves (ai = 0(e^)), is dictated by symmetry, 

dtA + sd^A = didlA + {ai + a2C)A-b\A\^A, (1) 

dtC = d2dlC - a^C + a^C^ + ih + h2C)\A\^ + h\A\^ + h^dMl^ + hdl\A\^ 
+hi{Ad,A - Ad, A) + h^id^iAd^A - Ad, A) + hsd,Ad,A. (2) 



The coefficients di,ai,a2,b in (1) are, in general, complex. The coefficients can be calcu- 
lated for a particular system through a perturbative, normal-mode expansion of the relevant 
variables with slowly varying amplitudes, 

$ = e{uA{x, t)e^(''*+^*) + c.c.) + e^{^rC{x, t) + ...) + h.o.t. (3) 

We consider a field C, that relaxes to a unique state in the absence of forcing by the wave and 
thus consider the contribution of the term proportional to 04, which could introduce a second 
branch through a transcritical bifurcation, to be of higher order. We furthermore focus on the 
affect of advection by the traveling waves on C by retaining only the gradient coupling in (0) 
and eliminating the dispersion in (|l]). In the case of thermal binary fluid convection it turns 
out that the term proportional to /ii, which is formally of lower order than the gradient term. 



only contributes at higher order |T6[. A gradient coupling of counterpropagating, traveling 



waves to a slowly varying field as a model for electroconvection in nematic liquid crystals 



has reproduced localized wave trains or 'worms' seen in experiment jT^. In systems with a 
true zero mode (conserved quantity) certain terms like those involving hi,h2, hs and /15 are 
not allowed and the gradient coupling is the relevant forcing |]T|, [§|. 

With the abovementioned simplifications, the system can be rewritten as, 

dtA + sd,A = dlA+{a + C)A-\A\^A, (4) 

dtC = 6dlC-aC + hd,\A\^. (5) 



Equations (|^, ^ describe the dynamics of a traveling wave without dispersion, which 
advects a real, slowly varying field C. The strength and sign of advection is given by the 
coefficient h. No homogeneous scaling is possible in these equations due to the presence 
of gradient terms which contribute to lower order in the limit of large-scale modulations. 
The scalings d, = 0(e), dt = O(e^) require s = 0(e), h = 0(e) for strict vahdity of (§,|) 
indicating small group velocity and weak coupling. Modulations on a longer length scale 
(dx = 0(e^)) allow for 0(1) values of s and h although the diffusive terms now contribute at 
higher order. The system thus becomes hyperbolic in this limit. Counterpropagating waves 
near onset have been studied in the hyperbolic limit ||1^, and can exhibit both sub- and 



supercritical secondary bifurcations as well as more complex dynamics. In this paper we 
choose not to reduce the equations (§,^ to a simpler form in a distinguished limit by fixing 
the scale. Instead we retain all terms as 0(1) quantities and consider the equations as a 
phenomenological model of travehng waves coupled to a critically damped mode. 



Our goal is to identify instabilities of waves in (^^ and characterize their nonlinear 
evolution. We note that in the absence of coupling {h=0) or in the limit of rapid decay of 
the C-field {a —>■ oo) (^^ reduce to a single Ginzburg-Landau equation which exhibits a 
long- wave phase instability for waves with wavenumber g^ > |. We shall refer to (|^,§) as 
the extended Ginzburg-Landau equations (EGLE). 

The complex amplitude A in (^ can be described by its real magnitude and phase 
A = Re^^. Thus, in a frame moving with the waves (§,|^) can be decomposed into the three 
following coupled equations, 



dtR = dlR + R[{a + C) - R^ - {d^e)% (6) 

R^dtO = d^[R^d,ei (7) 

dtC = 5dlC -aC + sd^C + 2hRd^R. (8) 

To determine instabilities of the traveling waves we linearize about the wave solution by 
considering an ansatz, 

R = y/a-q^ + e^e*P"+"*, (9) 

e = gx + e^e'P"+"*, (10) 

C = £^6*^^+"*. (11) 

We first consider instabilities in the longwave limit by expanding the growth rate a for 
small p. The three eigenvalues, corresponding to the real amplitude, phase and C-field 
respectively, are, 

a„„,(p) = -2R -2j^^^^^P 

[4i?V - 2shR^a + (R' + 2q') {2R^ - af] ^, ^ ^^^3^^ ^^2) 



i?2(2i?2_a)3 

a,,ase{p) = -DEp' + 2'^,p'-2^^^^^^p' + Oip% (13) 

i(2R^(s + h)-sa) 
ajMp) = -« + ^^^3^^ P 

l-AR'P + 2shR'a + i2R '-aY6] ^, ^ ^^^,^_ ^^^^ 



(2i?2 _ a) 
Here R^ = \/a — q^ is the amplitude of the plane wave and, 

^E = ^ — • (15) 

The parameter /? = h{s + h) is introduced for simplicity; it will prove to be useful in 
characterizing the type of instability. 



We notice that to leading order, the eigenvalues corresponding to the amplitude (|I^) and 
advected field (|1^ are negative for finite values of R^ and a. The critical mode, corresponding 
to the phase, exhibits a longwave instability when the diffusion coefficient De changes sign; 
this occurs at the Eckhaus curve g^ = |. However, the coupling of the advected field 
now raises the possibility of the quartic order term balancing the quadratic term in (12) 
as De -^ 0. This would indicate a small-, yet finite- wavenumber instability. Indeed, a 

2 2 

necessary condition for this phase instability to occur is that P + ^j^ < 0. Figure 1(a) shows 
the eigenvalue of the phase mode which exhibits this finite- wavenumber instability, i.e. 
sufficiently small wavenumbers are damped. In Figure 1(c) we see the corresponding linear 
stability diagram. Here plane waves become linearly unstable already before the Eckhaus 
curved is reached. However, we note that near onset of the traveling waves, i.e. as a — * for 
finite a, (^, ^ can be rescaled, and C adiabatically eliminated to yield a single Ginzburg- 
Landau equation. In this limit the contribution of the C-field is formally of higher order and 
we recover the Eckhaus instability (cf. Fig. 0(a), below). 

^"^ 2 2 

Next consider waves which are phase-stable in this regime {De > 0{p^) or (3 + ^^ > 
or both). In the limit as a — > 0, the eigenvalue corresponding to the C-field can become 
positive through a balance between the leading order term and that at 0{p'^). Consider, for 
example, the eigenvalues at bandcenter (q=0), 

a,„,(p) = -2. - 2-i^p - 1^°'^^ - ^;''"" + f - "'% ^ + Oij,% (16) 

[la — a) [2a — ay^ 

CTphaseip) = -P^ (17) 

, , i(2a(s + h) — sa) \—4a'^d + 2shaa + (2a — a)^6] o ^, qx,,^x 

cTf^eidip) = -^ + ^-^ S -P-- ^^ r ^^p' + Op3.18 

[2a — a) [2a — a) 

While for finite a as a — > 0, both the amplitude and phase eigenvalues are stable the 
field-eigenvalue can become unstable if —4/3 -t- 8aS < 0. This condition can only be satisfied 
for (3 > and implies additionally that there is an upper bound on the control parameter 
a for fixed f3 and 6 above which this instability cannot exist {abound = ^)- Such a case 
is seen in Figure |l](b) where the eigenvalue corresponding to modulations of the C-field 
indicates positive growth for finite-wavenumber perturbations. The linear stability diagram 
Figure |l](d) confirms the existence of the linear instability over a finite range of the control 
parameter a, resulting in a small stability-island near onset. 

As P passes from negative values through zero, the linear stability diagram given in Figure 
|l|(c) changes continuously to that shown in |l|(d). The intermediate regime is not captured 
analytically in the longwave expansions we have used thus far and must be investigated 
numerically. Figure ^(b) shows how the region of linear stability of the traveling waves 
changes as P passes through zero. The stability boundary, outside of which the waves 
are phase-unstable to finite-wavenumber modulations for /3 < 0, develops a bottle-necked 
region which pinches off as P becomes more positive, resulting in two separated plane-wave- 
stable regions (cf. Fig.|l](d)). This change is continuous, and there are values of P for 
which numerical simulations have revealed both phase- and amplitude-instabilities. That 
is, amplitude-modulated waves can arise whose underlying phase undergoes slips until a 
phase-stable wavenumber is achieved. 
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Figure 1: For all 4 figures h=l, a=0.02, 5=1.0 In (a) s=-1.5, q=0.425, a=1.0. This is a 
shortwave phase instability. The other two eigenvalues are large and negative. In (b) s=1.0, 
q=0, a=0.5. This is a shortwave instability of the C-field. The phase-eigenvalue is marginal 
for p=0 due to translation symmetry while the amplitude-eigenvalue is large and negative, 
(c) is the linear stability diagram for the values of the coefficients given in (a). The dotted 
line is the Eckhaus curve, given here for comparison, (d) is the linear stability diagram for 
the same coefficients as in (b). 



3 Phase Instability 

To facilitate the analysis of the phase instability, we make use of the large, negative eigenvalue 
of the magnitude of A to derive coupled phase and field equations in a longwave limit. If 
the longwave limit is taken at finite a both the magnitude and the advected field can be 
adiabatically eliminated and, to leading order, one obtains the usual phase equation with 
the usual Eckhaus stability limit. The effect of the C-field can be captured, however, in a 
distinguished longwave limit in which a — ^ as well. Then the amplitude can be adiabatically 
slaved to the dynamics of the phase and C-field, given that we are sufficiently far from the 
neutral stability boundary (i?^ ^ 1). 
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Figure 2: (a) Linear Stability diagram near threshold of the traveling waves (a ^ 1). The 
stability boundary approaches the Eckhaus curve as a — ;> (arrow), (b) Linear stability 
diagram with h=1.0, a=0.02, 6=1.0, s=-l.l, -0.6, -0.5, -0.452. In this intermediate regime 
both phase and amplitude instabilities can set in. 



We consider an ansatz of the form A = R(T,X)ei'^^'^''^\ where X = ex,T = e^t, and 
additionally a = e^a^- This scaling allows for large variations in the phase, although slow 
in time and space. To leading order we obtain an algebraic relationship relating the slaved 
amplitude to the C-field and local wavenumber. 
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{dx(t>f. (19) 

Combining equations at orders e and e^ yields the longwave equations for the phase and 
the C-field, 



(9t0 
drC 



DEdl<P + 



d 



'xyjux 



C 



5 die 



i?2 

tt2C+ (S 



h)dxC-hdx{dx(t>f. 



(20) 
(21) 



The phase-diffusion coefficient De is defined as before (cf. (0)) where E? is now given 
by (|19D and the local wavenumber is now given by the gradient of the phase, q = dx(p- As 
before, we linearize the longwave equations about the plane wave state with an ansatz. 



C 



qX + e0e^^^+'^^, 



(22) 
(23) 



The linearized longwave equations yield a complex, quadratic dispersion relation. As for 
the full EGLE we expand the growth rate a in powers of the modulation wavenumber p, 



0-phase{p) 
CTfieldip) 



-DeP^ + 2 



ihq^ 



-P 



-OL2 + i{s + Wjp - 6p^ 






/ + 0(p^ 



0{p^ 



(24) 
(25) 



As expected, phase modes with infinitesimal wavenumber become unstable when De 
changes sign. However, according to (0) for /? < phase modes with finite wavenumber p 
become unstable already for De > if 



p" > 2De 






(26) 



Thus, although we have identified the instability in the longwave limit, it is a shortwave 
instability and to determine the value of the critical modulation wavenumber p for which a 
first passes through zero we must solve the dispersion relation exactly. 

The action of the linearized operator of the longwave equations on (^, ^3D yields a 2x2 
matrix, the determinant of which must equal zero for a nontrivial solution to exist. If we 
consider the real and imaginary parts of the growth rate a = a-r + iuj, this dispersion relation 
can be written as, 



where 



a^ + ciar + C2 = 0, 
Csar + C4 = 0, 



(27) 
(28) 



Cl 
C2 
C3 

C4 

At criticality (a^ = 0), 



{pI{De + S) + 02) 



-LU 



up{s + h) + plDEiplS + 02) 



2lj — p{s + h), 
uj{p'^{DE + S) + a2) 






{DE{P-h') + h'). 



and 



^Hopf 



yield the Hopf-frequency, 

_ p3 {De{(3 - h^) + h^) 
~ h {p\De + 5) + ^2) ■ 



(29) 
(30) 
(31) 

(32) 



(33) 



From equations ( P7| ) and ( POJ ) we arrive at a bi-cubic polynomial in p, the roots of which 
give the critical modulation wavenumber for the oscillatory instability. 



-p\De{(3 - h') + hy+p'(5{DE{(3 - h') + h'){p\DE + 5) + a^) 
h^DEiSp^ + a2){p\DE + 5) + a^f = 0. 



(34) 



Solving for p directly in the sixth-order polynomial is not analytically feasible, but we 
note that the parameter (5 is merely quadratic in (31) and so seek a solution for /3, 
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1,2 



h\l~DE){p\6-DE) + a2) 

2DE{p^6 + a2) 
p{p^{DE + S) + a2 
'' 2DeP^{p^S + 02) 



/iV(l - DeY - ^DUp^ + aa) 



(35) 



It can be shown that cr^ > only if the real quantity (3 = h{s + h) is in the range f3i < 
(3 < (32- Thus, the finite- wavenumber instability can only arise if /?i_2 are real, requiring the 
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Figure 3: (a) Stability boundary as a function of the group velocity s. The boundary is 
bounded to the left by the stability condition derived from the dispersion relation of the 
longwave equations and to the right by the Eckhaus curve. Here h = 1, a = 0.02 b = 1.0 
and s = -1.5, -3.0, -5.0. (b) Dependence of the critical modulation wavenumber on the decay 
rate of C, a. The thick line is numerically calculated from the full dispersion relation while 
the thinner line is p=y^. ( ^OD is thus valid only for small enough a. (a=1.0, h=1.0, s=-1.5, 
5=1.0) 



discriminant to be positive, h?p'^{l — DeY > ^D\{p'^5 + a2)^- This can be interpreted as a 
constraint on the modulation wavenumber p. 
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(36) 



This implies the additional constraint h'^{l 



q'> 



DeY > 16-D|;(5a2 which can be simplified to 
a 



\h\ 



(37) 



2Va^ 



This somewhat intricate analysis of ( |5^ can now be summarized concisely. A sufficient con- 
dition for stability of the waves with respect to shortwave instabilities is that the underlying 



wavenumber of the waves be in the band given by q^ < 



3+- 



-, which is narrowed with 



respect to the traditional Eckhaus band. Thus, with respect to phase instabilities, the band 
of stable wavenumbers is bounded below by (E^) and above by the Eckhaus curve. 
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(38) 



If (|3^) is satisfied there exists a range in [3 given by ( p5D for which the plane wave is unstable. 
For a given j3 from this range the destabilizing modulation wavenumbers are given by the 
condition that the associated interval [/9i,/32] include that value of (3. 

In the limit of zero coupling {h — > 0) or large decay rate of the concentration field 
(q;2 — > C)o) the Eckhaus band is recovered in a consistent manner. A linear stability diagram 



of waves as obtained from the full EGLE (|I|, ^) is shown in Figure ^(a) for various values of 
the group velocity s. It can be seen that inequality (|38|) is satisfied. 



3.1 Weakly Nonlinear Analysis of Phase Instability 

To examine the weakly-nonlinear behavior of the phase-instability one must solve (|3^) for 
given values of the system parameters. In the general case this yields a numerical value for the 
critical modulation wavenumber p which can then be used in a weakly-nonlinear analysis. 
The coefficients of the resulting amplitude equation must then be evaluated numerically. 
However, if we restrict the wavenumber of the waves to lie on the curve representing the lower 
bound of existence of the shortwave instability, we can obtain analytical values for the critical 
modulation wavenumber and Hopf-frequency and consequently also for the coefficients of the 
amplitude equation. Along that curve. 
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It is good to keep in mind that these values fix a unique value of the parameter (3 = h{s+h) 
given by (^) and are in this sense restrictive. In addition, (|5^-^T]) are obtained in a longwave 
analysis, which requires that p, a <^ 1. Figure ^(b) indicates the values of a for which (|39| - ^TD 
are valid. 

We expand the phase and C- field in small-amplitude, normal modes, where e — - / ^°~°°^ 
is a measure of the distance from the bifurcation point (oq, Qcr), 
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(43) 



The complex amplitude ^40 of the unstable mode evolves on the superslow time scale 
r = e^T. Inserting this ansatz into (pO| , pT]) and solving order by order leads to a solvability 



condition at 0(e) which yields a differential equation for Aq, 

drAo = -'E{ao, a2, S, h)a2Ao + T{ao, a2, 6, h)\Ao\'^ Aq. 



(44) 



We are interested in the real part of the cubic coefficient, the sign of which gives the nature 
of the Hopf-bifurcation (forward/backward). The coefficient F is a complicated expression 
of 5, a2, a, and h, but we can expand it for small 0^2 to obtain an asymptotic approximation 
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Figure 4: Space-time diagram of the evolution of the magnitude \A\ leading to a phase slip 
(arrows), here indicated by the dip in amplitude. s=-1.5, h=l , a =0.02, 6 =1.0 



(a must be small from Figure H(b)), 
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Thus in the limit 0:2 — > 0, F > for all values of the coefficients, and the phase mode 
undergoes a backward Hopf-bifurcation. Solving the coefficient F numerically over a wide 
range of parameter values always yielded a positive real part. This analysis does not indicate 
whether or not the instability saturates at higher order or if the phase becomes undefined, 
signaling a phase slip. To investigate the nonlinear behavior of the phase, (|l|, ^ were 
integrated numerically using a linearized Crank-Nicholson scheme. It was found for the 
parameter values tested, that the waves outside the stability band underwent a phase slip 
when perturbed, relaxing to the plane-wave- stable band as seen in Figure ^, which shows a 
space-time diagram of the magnitude |^| of the waves. 



4 Amplitude Instability: Modulated Waves 

We now investigate the instability corresponding to the eigenvalue of the C-field passing 
through zero in (|T8|) which occurs only for /3 > 0. Figure ^(a) indicates that plane waves are 
stable at onset, becoming linearly unstable as the control parameter is increased until , for 
large enough values of the control parameter, they once again regain stability. In addition. 
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Figure 5: A linear stability diagram is given in (a) for h=l, s=l, a=0.02, 5=1.5. LW and SW 
denote longwave and shortwave instabilities respectively. In (b) the largest positive eigen- 
value of the system is plotted versus the modulation wavenumber at bandcenter. Shortwave 
instabilities occur as the control parameter a is increased from the stable region near onset 
of the plane waves and as a is decreased from the stable region far beyond threshold. 

the band of stable wavenumbers is bounded by the Eckhaus curve outside the range of the 
amplitude instability. Thus for /3 > two distinct instabilities are possible depending on the 
wavenumber of the waves and the distance from threshold. Figure |^(b) shows the relevant 
eigenvalue as a function of the modulation wavenumber for increasing values of the control 
parameter. The wavenumber of the fastest growing mode increases continuously over the 
range of values of the control parameter for which the waves are linearly unstable. 

At bandcenter, the system (|^, |^, ^) decouples in the case of unperturbed plane waves. 
Perturbing the plane-wave solution with the ansatz (j^, [1^, [Tl|) where q=0, introduces an 
order e^ coupling of the phase to the real amplitude through {d^OY, which does not appear 
in the linearized system. Thus it is sufficient to consider the following equations for the 
linear stability of the waves at bandcenter, 



dtR 
dtC 



dlR + R[{a + C)-R% 
bdlC -aC + sd^C + 2hRd^R. 



(47) 
(48) 



Linearizing about the plane-wave solution in the above equations yields, for spatially 
periodic modulations, a complex, quadratic polynomial in the growth rate a. As in section 
3, setting (T^ = gives a value for the Hopf-frequency, 

{2ao{s + h) + sp^) 

'^-^^=^(2ao + a) + (l + V)' ^''^ 

and a polynomial in p, the roots of which give the values of the critical modulation wavenum- 
ber, 



-p2(2ao(s + h) + sp^Y + sp^{2ao{s + h) + sp^)((2ao + «) + (! + 5)p^) 

+ (p2 + 2ao)(5p2 + a)((2ao + a) + (1 + 5^)2 = o. (50) 

If we are to solve (||),(0),(i]) to higher order, we must consider the couphng of the phase 
to the amphtude. We take an expansion of the form, 

Plugging this into (|^ yields the following series of equations for the phase, 

dA = d%, (52) 



^ n=l ' ' 



(53) 



Where n+j=k and k>0. Thus to leading order the phase satisfies the constant coefficient 
diffusion equation. Higher orders include a forcing term proportional to the gradient of all 
lower order terms which must all decay to zero for long times due to (|5^) . We conclude that 
the long-term behavior of small-amplitude instabilities at bandcenter will not be affected by 
the phase and consider (^, |4^) the relevant dynamics. 

4.1 Weakly Nonlinear Analysis 

In order to capture the dynamics of the instability near onset at both extremes of the 
unstable band in the control parameter, we carry out a weakly nonlinear analysis. Since 
solving the polynomial (^) for the critical modulation wavenumber is intractable we leave 
p as an implicit variable. 

We again introduce a slow time scale r = e^t and define the distance from onset e = 



y a ■ ^^ consider a small- amplitude expansion. 



C I 



^° ) Ao(T)e'(P-"+"-*) + c.c + h.o.t. 



and systematically arrive at a solvability condition at order 0(e'^) which yields a differential 
equation for Aq, 

drAQ = A(ao, a, S, h, s; p{ao, ..))a2Ao + n(ao, a, S, h, s; p{ao, ..)) | AqP^o- (54) 



The coefficients A and 11 in (p4D must be solved numerically by determining the value 
of p from ( pOD for given values of the parameters. For the values of the parameters tested, 
it was found that the secondary bifurcation encountered first as the control parameter a is 
increased from can be either forward or backward (cf. stability island including a=0 in 
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Figure 6: (a) Bifurcation diagram of modulated waves. Diamonds indicate minima and 
maxima of a single- wavelength modulation, squares and circles two and three wavelengths 
respectively. Solid symbols indicate the point at which the corresponding mode linearly 
destabilizes the traveling wave (triangles indicate four-wavelength instability although solu- 
tion branches are not shown). Parameter values are s=l, h=l, 6=1.5, a=0.02 with a system 
size L=125. (b) Solution with 2-wavelength modulation, (c) Switch of the bifurcation from 
backward to forward according to (0). (d) Comparison of the results from the amplitude 
equation (solid lines) to numerical simulations confirms this transition. Enough data points 
where taken for (A),(B) to indicate the presence of bistability and thus confirm the subcrit- 
ical nature of the bifurcations. s=l, h=l, a=0.02, and 6=1.1 (triangles), 1.2 (diamonds), 1.4 
(circles), 1.5 (squares) with a system size L=62.5. 



Fig. ^(a)). The bifurcation which occurs as a is decreased from the region of stability of the 
plane waves was, in all cases, supercritical. A bifurcation diagram for the modulated waves 
is given in Figure ^(a) for a system size of L=125. The particular structure of the cascade of 
bifurcations is a consequence of the finite system size, where each branch corresponds to a 
discrete number of wavelengths of the modulation. The bifurcation to modulated waves near 
onset of the traveling waves (for small a) is examined in detail in Figures §(c),(d) where the 
amphtude equation is compared to numerical simulation. As 6 is decreased for fixed values 
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Figure 7: Space-time diagram of \A\ showing the annihilation of one 'hump' by another 
leading to the formation of a single-wavelength modulated pattern from a two-wavelength 
pattern. s=l, h=l, a=0.02, 6=1.5, a=0.04 



of the other parameters, the bifurcation switches from forward to backward. Figure ^(d) 
confirms the validity of the amplitude equation. 

The dynamics of the modulated wave pattern is not trivial. The destruction or creation of 
a modulation wavelength during the transient dynamics before the stable pattern is achieved 
occurs due to particle-like interactions of the 'humps' of the modulations as seen in Figure 
^. The bifurcation diagram in Figure ^(a) indicates the existence of multiple branches for a 
single value of the control parameter. If the initial bifurcation is to a lower branch (say with 
3 wavelengths), two of the three 'humps' will collide after an initial periodically modulated 
state. The remaining two may again collide to form a single-wavelength modulation which 
is then stable. 

As 6 is decreased further, the subcritical branch of modulated waves, created in a saddle- 
node bifurcation as the control parameter is increased, extends further towards the onset 
of the traveling waves, eventually becoming bistable with the basic, conductive state (Fig. 
^(a)). As in the supercritical case, the modulated- wave branches arise from an instability of 
the plane-waves to finite wavenumber modulations. Modulation wavenumbers p which lie in 
the interval pi < p < Pu, where the bounds depend on the value of the control parameter 
a, will linearly destabilize the plane-waves. For a finite system of length L, there will exist 
modulated-wave branches for all n such that p = ^^ lies in this interval. As n is decreased 
the solutions increasingly take on the characteristics of localized objects (Fig. ^(b)). In 
fact, the distance between them can become arbitrarily large. Consequently, this distance 
can be much larger than the maximal wavelength — of modulations that linearly destabilize 
the traveling waves in a given system. Thus, these solutions do not bifurcate directly from 
the plane-wave branch but rather most likely arise in a tertiary bifurcation off a modulated- 
wave solution. For example. Figure § shows a series of modulated-wave branches for a system 
(L=250) where a single modulation (n=l) does not linearly destabilize the waves. Still, the 
single-pulse solution exists (dash-dot branch). 
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Figure 8: (a) Bifurcation diagram of pulses; the subcritical branch extends over the con- 
ductive state. The solid symbols indicate where the modes of the corresponding branches 
become linearly unstable, (b) Profile of the pulse solution with h=l, s=l, a = 0.02, d = 0.5, 
a=-0.1. The solid line is the amplitude \A\ and the dashed line the C-field profile. 



5 Conclusion 

Systems similar to (3), (4) have been derived in various contexts, including a model of trav- 
eling interfacial waves in two-layer Couette flow [^ where a = due to conservation of 
the additional zero-mode, given by the position of the interface. This system was shown to 
exhibit a phase- instability leading to a phase-slip in agreement with experiment ||l|. A 
small-amplitude model for traveling waves in thermal, binary-mixture convection consists of 
complex Ginzburg-Landau equations for counterpropagating waves coupled to a slowly de- 
caying mode representing large-scale variations of the concentration field @]. In this system 
the traveling waves arise in a backward bifurcation. Localized pulse solutions arising from 
instabilities of this system were identified and many of their properties characterized ||16| . 

In this paper we have shown that the advected mode can cause the phase-instability 
to occur at finite wavelength and can introduce an additional amplitude-instability. The 
latter is identified as the origin of pulse-solutions seen in earlier works. They arise when 
the secondary bifurcation to modulated waves is sufficiently subcritical to lead to bistability 
between the modulated wave and the basic state. A similar bifurcation may explain the 
appearance of localized "worms" in electroconvection in nematic liquid crystals, seen again 
before the primary instability to supercritical traveling waves. Indeed, a model of traveling 
waves coupled to a slowly decaying field has been invoked to describe the worm dynamics, 
attaining qualitative agreement [14 . 

Finally, it would be worthwhile to consider the effect of dispersion in (^) on the dynamics 
discussed in this paper. In some cases dispersion has been shown to have a significant effect 
on the linear-stability properties of waves coupled to a additional mode. A model similar to 
(|, H) used by Bartelet and Charru [|l|, with a = and added dispersion, seems to explain 
qualitatively the linear-stability properties of waves observed in experiment. The inclusion 
of dispersion also allows for a determination of the effect of the additional mode on the 
Benjamin- Feir instability. In some cases this effect is important as, for example, in the 4- 
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Figure 9: Bifurcation Diagram with /i=1.0, s=1.0, a=0.02, (^=0.5, g=0 and L=250. Each 
branch represents a modulated-wave solution with the number of modulations given. As this 
number n decreases, the modulated- waves become more like localized traveling pulses. 

species Oregonator model [Q of the BZ reaction where the Benjamin- Feir instability can be 
preceded by a finite-wavenumber instability. 
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6 Appendix A: Derivation of Amplitude Equation for 
Phase Instability 



Here we outline the derivation of the amplitude equation ( PI ) which describes the slow- 
time evolution of modulations of plane waves in (|20| , ^I]) . We consider finite-wavenumber 
modulations pcr of waves with wavenumber Qcr as given by ( ^OD and (p9|). We first expand 
the phase and C-field as, 



C 



QcrX 




e$o + e^^i + • • • , 



(55) 



where e 



2—^ and ao is the value of the control parameter at onset of the fastest-growing 



mode. We rewrite (|20| , |2l|) symbolically as 



a ^ )+Ar(0,c) = o, 



and expand the linear and nonlinear operators in e. 



(56) 



(57) 



Af = e^N, + e^N^ + ..., (58) 

De = D + e^D2, (59) 



with 



"^^ V -2%.5i <5ai - a^ - a2 + (s + /i)ax J ' ^^"^ 

Li = 0, (61) 

iVi = {KiCd\(t) + K2dx(t)d\(t) + K^CdxC + Kidx(t)dxC, -2hdx(j)dj^(t) f , (63) 

' V K^CdxCdxc/^ + K.oCdxC, ii J 

(64) 



KgCdxCdx<p + KioC'^xC, 

where the coefficients are given by, 

a2qc 



R = \/ao-ql.,R2 = — ^, (65) 

O = ^^■D, = ^fi-D), (66) 

K. = ^^.A'. = ^P-3), (67) 

A'r = ^(d + «o), /^8 = |p(3 + ^), (70) 

The slow time is given by r = e^T. At order e we recover the homogeneous, linearized 
system, 

I^o^o = 0. (72) 

The zero-mode $o we take, as in 3.1, to be of the form, 

$0 = [ ^° ") ^p(7-)e^(p-^+--^) + c.c. (73) 

The eigenvector {^q,Co)^ is determined by solving ([7^), 



$0 \^ /^ ^^crff 

Co y V ^Pcr + ^^cr 



2 ,--■,, ) (74) 



The left-null eigenvector is found by solving the adjoint problem, 



At Ofe^) we have 






Lo$i + iVi($o)=0. (76) 



The nonlinear term A'^i contributes both spatio-temporally homogeneous modes and 
modes of twice the critical modulation. We thus take $1 to be of the form, 

^^=f ^12 \ ^2g2.(p..X+...T) + c.c + (" ^^^^° ] \Ao\'. (77) 

With this ansatz, (^) reduces to two linear, algebraic systems. The system corresponding 
to the forced second harmonic is 

4Dpl - 2iu,r 2zper|f W $12 

8hqcrPlr -^hlr " " + 2«(s + h) - 2iuJcr J V C'la 

Ki^oCopi + iK2^lpi - iKsC^pcr + K^^oCopl^ 
~2ih^lp% 

and that for the homogeneous forcing 



0, (7^ 



-1 \ / $10 \ / -K,pU%Co + $oQ) + K,pI{%C, + <l>oCo*) 
-a j I Cio J ^ V 



0. (79) 



Note that the homogeneous forcing causes a response in the phase which grows linearly with 
time, i.e. it changes the frequency of the wave (cf. (ffTD). 
At O(e^) we have 

Lo$2 + L2^o + Ni($o, $1) + Ni($i, $0) + N2($o, $0, *o) = 0. (80) 

At this order the critical Fourier mode appears through the nonlinear forcing terms and 
the action of the linear operator L2 on the eigenvalue of the linearized operator of the original 
system. In order for a solution to exist, these contributions must lie orthogonal to the left 
nuUspace of the operator Lq. We thus obtain the solvability condition, 

($Jj*, L2$o + Ni($o, $1) + Ni($i, $0) + N2(#o)) = 0, (81) 

where the inner product is defined as an integral over one period of the critical Fourier mode 
and the scalar product of the scalar vectors in $, C-space. The solution of ( pT]) results in a 



differential equation (44) for the complex amplitude Aq, the expression for which is extremely 
long and involved. 



7 Appendix B: Derivation of Amplitude Equation at 
Bandcenter for Amplitude-Driven Instability 

Here we follow the same procedure as in Appendix A in outlining the derivation of the Landau 
equation for the amplitude-driven instability. The calculation is carried out at bandcenter, 
i.e. q = 0. We perturb the basic, plane-wave state in equations (|i7| , ^ ) by expanding the 

real amplitude R and real field C in e = \ h^f^ where a^ is defined as in Appendix A, 

^) = (^)+.*o + e^*i + .... (82) 



We rewrite (E^ HSF) symbolically as. 



^{c^ +Ar(i?,C) = 0. 



We now expand the linear and nonlinear operators in e. 



with 



' h^d, -dr 



Co J ~ \plr + i^cr + 2ao 



(83) 



£ = Lo + eLi + e^L2 + ..., (84) 

AT = e^Ni + e^N2 + ..., (85) 



J _ I di-dt-2ao ^ \ , . 

"^0 ^ 2h^^d., 5dl-dt-a + sd., ) ' ^""""^ 

Li = 0, (87) 

-2a2-dr ^^ 



A^^ = [RC-?,^R^, 2hRd^Rf, (89) 

N2 = ( -i?^ 0)^. 

(90) 
At order e we recover the homogeneous, linearized system, 

1^0*0 = 0. (91) 

The zero-mode ^o "we take, as in 4.1, to be of the form, 

^0= ( ^^ ^g(^)e^(p-^+--*) + c.c. (92) 

The eigenvector (i?o, C'o)"^ is determined by solving (|9TD , 

(93) 



The left-null eigenvector is found by solving the adjoint problem, 



\ Per + 2^0 — ioJcr J 



At Ofe^) we have 



Lo*i + A^i(*o)=0. (95) 

The nonlinear term Ni contributes both spatio-temporally homogeneous modes and 
modes of twice the critical modulation. We thus take ^i to be of the form 

^^=f ^|2 ^ ^2g2.(p..X+...T) + c.C + r ^1° ] \Ao\'. (96) 

With this ansatz, (^) reduces to two linear, algebraic systems. The system corresponding 
to the forced second harmonic is 

-Apl. - 2iucr - 2ao ya^ A / R12 

Ahy/oQiPcr —4:dpl^. — 2iuJcr — a + 2ispcr J \ C12 



2ihpcrRo 







0, (97) 



and that for the homogeneous forcing 

-2ao y/a^ \ f RiQ \ f RqCq + RqCq - 6^/a^\Ro\'^ \ _ n fnQ\ 

Solving these systems yields the following expression for the scalar quantities, 

Ri2 = ^, (99) 

"12 

a. = ^, (100) 

"12 
"^12 = - - A/a^(4ao(5Pcr + 4<^Pcr + «0« + Plr<^ + '^■UJcrSPcr - 2cJ^^ + Z [ASpl^UJcr 

+aujcr - 2aospcr - 2spl^ + 2aoUcr + 2UcrPlr + 2ao/ipcr]), (101) 

ni2 = -Aaohpcr ( -uJcr + i[aQ + -plA ) > (102) 

di2 = -86pl^ - 2{a + 2ao6)pl^ - a^a + 2uJcr{uJcr - spcr) 

i[4:spl^ - 4u;„-(l + S)pI^ + 2aoPcr(s + h) - Ucri2ao + a)], (103) 

Rio = -V«o + ^(Pcr + 2ao), (104) 

Cio = 0. (105) 

At cubic order projection onto the left null eigenvector yields the amplitude equation for 

drAo = I {pcr ) {plr + iuJcr)Ao + l {pcr ) aAo 

[Ci2V^ - Sal^^ + 2R^2{pI - ao) + 2R,o{pI - ao + ^u;cr)]|^^^, (106) 
{{pI. + 2ao)^ - cug^ - 2i[ucr{plr + 2ao) + happcr ]} 
[{pI,. + 2ao)2 - a;2J2 _^ 4[u;^^(p2^ + 2ao) + /laoPcrj 



^(pc.) = 2^hp^r ])7 :::z :\. , :r ^ t r; t tz'X • (lo^) 
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